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1. Introduction. 

The uncertainty of the primary cosmic ray composition 
at 10' 4 -10 16 eV is well known to make the study of the 
nuclear interaction mechanism more difficult. Experimentally 
considering, if one can identify effectively the family 
events which are produced in low heights, then an event 
sample induced by primary protons might be able to be 
separated. It is undoubtedly very meaningful. In this 
paper we simulate the family events under the condition of 
mountain emulsion chamber experiments with a reasonable 
model. The aim is to search for the dependence of some 
experimentally observable quantities to the interaction 
height. 

2. Method. 

The model used is: proton incidence, total inelastic 
cross section rising with energies ( O =* E°‘° 6 ), Feynman 
scaling in fragmentation region holding and mean transverse 
momentum increasing with energies ( <Pt> E )• The 
electromagnetic cascade and the multiple Coulomb scattering 
are also treated by Monte— Carlo method. Because the stress 
point is on the study of the qualitative charaters of the 
production height dependence on the transverse distribution 
of family events, the model used may be a suitable one. 

The observation height is assumed at Mt. Kanbala (520 
g/cm^). The incident axis is taken uniformly on 40x50 cm 2 
X-ray film, and the particle outside this area are ignored. 
Hadron showers are recorded by the same method as experiment 
with chamber thickness of 28 c.u. Pb. The criteria of a 
family are 2(E y +E£”) (=2E)>30 TeV, n y ^ 4 and n h ^ 0 (all 
symbols are in the common meaning). 546 family events with 
IE £100 TeV are analysed. 

The (R) and <ER> of a group of testing events with 
<Pt>=330 MeV/c are shown in Table 1 to compare with 
M.Shibata's results (in parentheses) using model PRS*'* # The 
difference between them may be reasonable owing to the 
different observation levels. For families with IE £ 100 TeV 
(Ewiia =4 TeV), we have <R>=2.75 cm which is close to the 
experimental value 2.5 cin 01 • 
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Table 1 <R> 

and <ER> 

of family 

events 

E 

(TeV) 

30-50 

50-100 

100-200 

200-500 

<R> 

(cm) 

3.20 

3.45 

2.82 

2.43 

(3.3) 

(3.1) 

(2.5) 

(2.3) 

<ER) 

► (TeV. 

cm) 12.5 

14.1 

12.2 

10.7 

(13.6) 

(13.7) 

(11.7) 

(10.6) 


3. Average interaction times, mean interaction height 
and purity. 

The average interaction times, l\!;nt , for all observed 
particles of each family event is calculated. In Fig.l is 
shou/n the N ^rvt distribution with mean value 3.74. The 



Fig.l 


Fig. 2 



relation between IMint and the first interaction height, H 1$ t 
is seen in Fig. 2. It is obvious that the events having both 
single interaction only and low interaction height are very 
rare. 

For each event with multi-time interactions, imagine 
that there exists a main interaction. We devide the 520 
g/cm 2 atmosphere into 21 layers with step 25 g/cm 2 * Define 
the layer in which the produced particles contributes the 
most fraction of family energy 2 E as the main interaction 
layer and the corresponding height as main interaction 
height (Hmaln). The fraction mentioned above is defined as 
purity of the family (abbreviate as p hereafter). 


The distribution of main interaction height is rather 
flat, slight favorable to lower height (or larger Hmaln) 

(see Fig. 3). It is noticed that the events with larger 
H main are not always very pure (Fig. 4). The low-height event 
which we are interested in, of caurse, should those with 
lower main interaction heights and high purities. 


R and ER. 

The analysis shows that <R> and <ER> of family events 
are significantly dependent on Hist (Fig.5): events with 
larger <R> and <ER> have smaller Hist , events with larger 
H^st have smaller <R> and ^ER>, but incorrect inversely. As 
to the dependence of TR> and <ER> on H truxXrx f it is insens i- 
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tiv/e (Fig, 6). 

The rejuvenation treatment , or taking R as the dis- 
tance to event axis, or taking <R> as the energy weighted 
mean R, all induce the <R) and <ER> with the same qualitative 
characters as shown in Fig, 5 and 6, 

"Decascading" treatment with full efficiency gives weak 
dependence of <R> on H ma;n (fig. 7). Besides, ^ER>, mean 
energy (<Ey> ) and maximum energy (Ey?** ) of y-rays of a family 
event show the same tendency as <R>. But at the other hand 
y-ray number of a family, ny , shows opposite tendency. In 




Fig, 8 is shown the dependence of (Ey> on H Wft ; n • When we 
analyse the relation between the composite quantity Cs^R>-{Ey) 
Ey* x /ny 1 and H waJn (Fig, 9), it is seen that some low height 
events are condensed near C=0. Taking a C value cut (say 
C=5), the event sample selected will be dominant for low- 
height events. 
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Fig. 7 


Fig. 8 


For the families not treated by decascading procedure, 
the composite quantity A ~ ny /( <R>-<ER>-<Ey>-Ey fty ) can also be 
used to separate the the concerned events . J For example, 
there are 28 events satlsfing the conditions A<0.05 and <TR> 

< 3 cm, among which 12 have Hw#;n ^ 450 g/cm , 7 events have 
both Hw,a; n > 450 g/cm and p>0.7. Considering the whole events 



satisfing these two condi 
tions are not very rich (Fig. 
4), this method may be use- 
ful. 

The condition <ER><10 
TeV.cm and <Ey/>12 TeV can 
also separate about one half 
events having H mo;n > 450 g/cm 
and p>0.8. 


Hmain (9/cmh 

Fig. 9 


Summary. 

Due mainly to the multi- 
interactions in thick atmos- 


phere target above the mountain emulsion chamber, the depen- 
dence of the observable quantities of family events on the 
main interaction height Hm*:n is strongly smearing. However, 
weak dependence on H w< «n still exists in <R> , <ER> , <Ey> , n> 
and Ey** etc. A small sample with low— height and high purity 
is not impossible to be separated by suitable composing and 
cutting of these quantities. 
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